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Abstract 

We propose a scenario where Dark Matter (DM) annihilates into an intermediate state which 
travels a distance A = v/T on the order of galactic scales and then decays to Standard Model 
(SM) particles. The long lifetime disperses the production zone of the SM particles away from 
the galactic center and hence, relaxes constraints from gamma ray observations on canonical an- 
nihilation scenarios. We utilize this set up to explain the electron and positron excesses observed 
recently by PAMELA, ATIC and FERMI. While an explanation in terms of usual DM annihilations 
seems to conflict with gamma ray observations, we show that within the proposed scenario, the 
PAMELA/ ATIC /FERMI results are consistent with the gamma ray data. The distinction from 
decay scenarios is discsussed and we comment on the prospects for DM production at LHC. The 
typical decay length A 10 kpc of the intermediate state can have its origin from a dimension 
six operator suppressed by a scale A ~ 10 13 GeV, which is roughly the seesaw scale for neutrino 
masses. 
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I. INTRODUCTION 



An exciting and very plausible possibility is that dark matter (DM) interacts non- 
gravitationally. A very active ongoing program exists to search for such interactions both 
directly and indirectly. Indirect evidence could arise from cosmic ray signals originating from 
areas of high dark matter density — the galactic halos and cores. Decaying or annihilating 
DM can act as an additional source of cosmic ray fluxes on Earth, and could be seen as an 
excess above the expected flux. The channel in which the excess could show up depends on 
how DM couples to standard model (SM) particles. 

DM decays and annihilations into charged particles inevitably lead to enhanced gamma 
ray fluxes due to internal and final state bremsstrahlung. In both cases the gamma ray 
spectrum will be correlated with the energy spectrum of charged particle, and will be the 
same (up to overall energy scale) for the same final states. However, the angular distribution 
of gamma rays in the sky will differ between the two cases. Since the annihilation rate scales 
with p 2 (p being the DM density), the gamma ray signal will be more highly peaked toward 
the galactic center due to the rise in the density profile for small r. The decay rate, on the 
other hand, scales linearly with p and therefore the gamma ray signal is less peaked. 

In this paper we consider an alternative scenario, where the DM annihilates into long 
lived particles (LLP). This setup will lead to a gamma ray signature which can interpolate 
between the signatures of DM decays and annihilation. The prolonged decay lifetime of LLP 
effectively smears out the distribution of gamma rays. To be concrete we will consider this 
scenario in the context of DM explanation of the positron excess found by PAMELA [lj], 
ATIC [2| and FERMI [jjj]. The PAMELA experiment which sees a rise in the positron-to- 
electron ratio in cosmic rays for energies between 10 and 100 GeV corroborates and extends 
the range of earlier data [4| , while ATIC sees a rise and then a striking fall near 800 GeV in 
the power-law for the combined electron and positron flux. The rise at lower energies, below 
the bump, corroborates previous data ^ while the bump at 800 GeV had never been seen 
before. Recent data from the FERMI satellite do not confirm this bump, while an excess 
of the electron and positron flux over the expectations from common cosmic ray models is 
still there. FERMI data indicate a smooth spectrum falling approximately as E~ 3 till about 
800 GeV. Above this energy data from the HESS telescope jy, 0] indicate a much steeper 
slope of the electron/positron flux. Although it seems quite reasonable to believe that this 
data has an explanation in terms of nearby young pulsars j§, 0, QiJ 11, 12, 13, 14, 15] (or 



some of the data could be spurious), it will serve as a useful tool in highlighting the features 
of the LLP scenario. In particular, the data exemplify a case where there is a need for an 
increased annihilation cross section (<jv) ~ 10~ 22 cm 3 /s (relative to the typical cross section 
(av) ~ 3 x 10 _26 cm 3 /s needed for Qbm ~ 0.2), while there is no evidence of excess in the 
gamma ray flux in the direction of the galactic center. Standard annihilation scenarios with 



such large cross sections seem to conflict 16|, [171, [l8|, [l9j with the results from the HESS 



telescope [20l . |21| for commonly assumed DM profiles. We will show that, by allowing the 
annihilation products to be long lived, the constraints from gamma rays can be avoided. 

To be specific, let us consider the situation where a DM particle x annihilates into an 
intermediate state </> which subsequently decays into standard model particles 

XX ^00^2SM2SM. (1) 
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As an example we will assume in the following that decays into muons, 



(2) 



though our arguments do not rely on the specific decay mode of 0. We concentrate on this 
mode since it is the relevant one for the PAMELA/ATIC data (while not changing our main 
conclusions other annihilation modes may be preferred when including new FERMI data 
0)- In contrast to @, 0, 0, H 0, |28| . [29( we now assume that is a long-lived particle 
(LLP), such that it will propagate over galactic distances. 1 For the moment we will not 
concern ourselves with whether or not the DM particle \ is a thermal relic. We will return 
to this issue at the end of the paper in section HVl In Sec. [IT] we present a general discussion 
of the LLP scenario and discuss the effects on the production of Standard Model particles 
from DM annihilations focusing on gamma rays. In Sec. IHII we apply this idea to the recent 
electron and positron flux data from PAMELA, ATIC, FERMI and HESS and show that 
they can be made consistent with gamma ray observations from HESS even for rather cuspy 
DM profiles. In Sec. |V]we summarize our results, discuss possible distinguishing features of 
the scenario, mention prospects to see DM at LHC, and speculate on a possible connection 
between the high scale physics responsible for the decay of and the suppression of neutrino 
masses via the seesaw mechanism. 



II. GAMMA RAYS FROM DM ANNIHILATIONS VIA A LONG-LIVED STATE 
A. General discussion 

Let a be produced at f = 0. Then the probability that it decays in a volume element 
d 3 r around r is 

- x —i>r, with (3, 

where T is the decay rate of 0, A the corresponding decay length in the laboratory frame, 
while the velocity of in the laboratory frame, f3 = v/c, and the relativistic gamma factor 
are 



m A/> V 7 V m x 

where m x ^ are the masses of \ and 0. In the last equalities above we used the fact that \ 
is non-relativistic. 

Let us now calculate the photon flux resulting from final state radiation of the process 
shown in Eqs. [1] and [2j The differential photon flux $ 7 in a solid angle AQ is given by 

^ - 2 <„> / *> f ds f^imj-^i -i . d j-< . . ( 5 ) 
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dE^ J An Ao.s. J 2 ?77^ 4n\ \f— r '\ 2 2ndE 1 dcos8 



1 Phenomenologically, annihilations into rapidly decaying particles are challenging to distinguish from pure 



annihilation 
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301 ] and thus we will not consider them distinct in this paper. 
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Here p(r) is the energy density of DM particles Xi ( av ) is the averaged annihilation cross 
section of \i an d we have assumed that \ is a Majorana particle, a choice we shall adopt 
throughout this work. The factor of 1/2 accounts for the fact that the incoming particles 
are identical, as it has not been included in the definition of (cxt>). d 2 N 1 /(dE 1 d cos 9) is 
the double differential gamma spectrum per decay. It can be obtained by boosting the 
isotropic spectrum in the rest frame (RF), where 



d 2 N^ 



dE^d cos 9 



ldN^ 



RF ^ dEry 



(6) 

RF 



The factor 2 in Eq. [5] accounts for the fact that we will have two decays per annihilation. 

In Eq. we use coordinates with the origin at the galactic center, and XX annihilate at 
the point r' and decays at r. Then 9 is the angle between the line of sight (l.o.s.) and the 
direction of motion of 0: 

cose = -t£zp- (7 ) 

where s is the unit vector along the line of sight. In the following we are going to compare 
the predicted photon flux to data from the HESS telescope, using observations of the galactic 
center (GC) [20| within a cone of solid angle 10 -5 . In this case we have 

p p2w pAxp 

A^gc = 2tt(1 -cosAVO = 1(T 5 , / dVt = / dip sin^#, (8) 

J An Jo Jo 



s 2 — 2r s cos ip , s = (cos-0, cosipsmip, sin 99 sin -0) , (9) 



where r Q ~ 8.5 kpc is the distance of the solar system from the galactic center. Furthermore, 
we use the HESS observation of the galactic ridge (GR) j2l|. This is a rectangular region of 
the sky with galactic latitude and longitude of |/| < Al = 0.8°, |6| < Ab = 0.3°, where the 
inner region corresponding to the GC defined above has been subtracted. For the GR we 
have 

i GR m AAlAb — Af2cc ~ 2.8 x 10~ 4 , / dtt w 4 / dl db , (10) 

J An Jo Jo 

I / l 2 + b 2 \ 

r| + s 2 — 2r Q s cosZ cos 6 , s « ( 1 — , l,bj . (11) 



To compare with HESS observations [21(, one additionally needs to subtract from the photon 
flux observed in the GR region |/| < 0.8°, |6| < 0.3° the average photon flux in the region 



0.8° < |6| < 1.5°, a procedure that was used for the background subtraction in [21 



Before we specialize to the case of non-relativistic <fi and proceed with our discussion we 
mention briefly some limiting cases of Eq. [5j First, we consider the limit r — ► 00 (or A — ► 0). 
In this case the integrand in Eq. [5] will be non-zero only for r = r'. Hence we can replace 
p(r') — > p(r), which can be pulled out of the d 3 r' integral. Using for the remaining d 3 r ; 
integral 

"dry 1 ^ _dN, 

2iTdE y dcos9 dE y } y ' 



4 



one obtains 



d$>. 



(av) c?iV 7 
47r dE^ 



dn I ds^-. (13) 

T^oo "^7 J AH Acs. m x 

As expected we recover the standard expression for the bremsstrahlung photon flux from 
annihilation XX ~^ + wy, apart from a factor 2, since in our case there are twice as 
many photons because of two decays per x annihilation (see Eq. [2]). 

Second, we consider the case of highly relativistic 0. Then one has 7 > 1 and (3 ~ 1, 
while the photons are collimated with within an opening angle of size I/7. This implies 
that there is only a contribution to the integral when (r — r') is aligned with the line of 
sight, and the d 3 r' integration reduces to another integration along the line of sight: 



(av) dN^ 



rclat. 



4tt dE. 



dtt 



■yJAn 



ds' / ds 



p\r>) 

771?, 



(av) dN 1 



4vr dE. 



dVL 



-yJAQ 



ds' 



l.o. s. 



771?, 



r (s'-s)/X 

A 

-s'/X 



(14) 



As a reference, we also give the expression for the photon flux from decaying dark matter 
(see, e.g. HU) 



dE^ 



r dAL 



decay 



4tt dE. 



dVL 



7 "'An 



p(r) 



ds 

l.o.s. m 



(15) 



where V is the DM decay width. The important difference between decaying and annihilating 
DM is that the DM density enters linearly in the first case, and quadratically in the second, 
which has important observational implications discussed below. 



B. The effective DM profile for non-relativistic intermediate state 



In order to simplify the calculations we specialize from now on to the case of non- 
relativistic 0, as the generic effects will not change for the case of a relativistic <fi. The 
effect of an intermediate LLP with a suitable decay length A will always be to flatten out 
the distribution of SM particle production with respect to the DM distribution, indepen- 
dently of its 7. This follows from the structure of Eq. [TH or the general expression Eq. |5l 
The general case would deserve a dedicated quantitative study which is beyond the scope of 
the present work. 

For non-relativistic one has 7 « 1, (3 <C 1, and the decay of <fi is isotropic. In particular, 
the photon spectrum is given by Eq. The spectrum can be pulled out of the integration 
and the effect is a smearing of the density distribution on the scale A. One may define an 
effective density distribution: 



Peff( r ) 



d 3 r'p 2 (r') 



I e -|?-r'|/A 

47rA \r — r'l 2 



1 f°° , , , o, „ L. ( r + r' 



2Xr 



dr'r' p 2 {r') 



Ei 



with the exponential integral 



Ei(x) = -PV 



00 e - t 



Ei 



dt- 



(16) 



(17) 
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where PV denotes principal value. Using the effective DM profile in Eq. [16] we can write the 
photon flux from a given solid angle AQ as 



d$ 7 _ (av) r p dN 7 
where the dimensionless J-factor is defined by 



JAfl (18) 



L[ d nf 

AO J AQ J Lo . s . r p| 



J = _i_ / d n I — (19) 



For the numerical calculations in this work we use p = 0.3GeVcm 3 and r = 8.5 kpc 



and for the DM density p(r) we always assume a NFW profile [32 



r (\ +r Q /r s 



2 



p(r)=Po^ > ( 2 °) 

r \ 1 + r/r s J 

with r s = 20 kpc. Eq. [18] has the same form as the photon flux from standard annihilations 
XX ~ * P + P~ + nj, apart from a factor 2 from the decay of two 0's, as mentioned above. 

In the non-relativistic case the two parameters and T appear only in the particular 
combination A = 7/3/T, see Eqs. [3]and[H Therefore, apart from the two parameters (erf) and 
m x of the standard annihilation scenario, we have now effectively one additional parameter 
corresponding to the decay length of in the rest frame of the galaxy. Numerically one has 



r r * \o.oi J vio 12 s 

with 5 = (m x — m ( j } )/m x 1 and r = 1/T. 

In the LLP scenario the source of SM particles from DM annihilations is proportional to 
Peff( r )> "while in the case of standard annihilations it is proportional to p 2 {r). This means that 
via the propagation of the intermediate state (ft over galactic distances we decouple to some 
extent the production of SM particles and associated gammas from the DM distribution. 
Fig. [1] shows the effective DM profile for various values of A. From the figure we find that 
for r < A we suppress p e &{r) with respect to p{r) since the had no time to decay yet, while 
for r ^ A we obtain a p e fr slightly larger than the original DM profile. The relative size of 
this over-production increases with A. Finally, at large distances r ^> A all 0's have decayed 
and we recover the NFW profile. 

One can understand this behavior qualitatively from the definition of p e ff in Eq. [16] From 
the first line of this equation it follows that for r ^> A the exponential is non-zero only for 
r ~ r' and one obtains p c fr(^) — > p( r )- O n the other hand, for r < A and for profiles 
p(r) oc r~ 7 (7 > 1) one has roughly 

2 , f d 3 r' f dndr' 1 

PeffW| r<<A oc J J oc (22) 

where the last relation follows just from dimensional analysis. Hence, the slope of the photon 
source term at r A is reduced by one power with respect to p 2 {r). Since for a NFW profile 
7 = 1, we find for r <C A that pl s (r) oc PNFw( r ) oc 1/r. Note that for DM decay the source of 
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r [kpc] 




FIG. 2: Interpolation between the SM particle source terms for DM annihilation and decay. The dashed 
curves show PnfwM// © f° r annihilation and /?nfw (j)// 9 © f° r decay, whereas the solid curve corresponds 
to Poff( r )//°0 f° r DM annihilation into a long lived intermediate state with a decay length of A = 10 kpc. 



gamma rays is proportional to p(r), in contrast to the p 2 (r) for annihilations. Therefore, for a 
NFW profile our scenario exactly interpolates between DM decay at r A and annihilation 
at r ^> A. This behavior is shown in Fig. [21 Note that the parameter A controls the absolute 
value of plfi at small r. For A ~ 40 kpc one would match exactly the profile for DM decay in 
the particular (somewhat arbitrary) normalization used in Fig. [2J However, for such large A 
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FIG. 3: J-factors for GC and GR (background subtracted as in 2l|) as a function of A. Dashed lines 
correspond to the J-factors for A = 0. 



A [kpc] 





1 


10 


100 


GC 


14200 


781 


101 


10.4 


GR 


1430 


174 


18 


1.9 



TABLE I: J-factors for GC and GR (background subtracted as in 2l|) for various values of A 



the asymptotic p 2 (r) behavior at large r is reached only at distances much larger than the 
size of our galaxy. We emphasize that the slope at small r is independent of A, and hence 
it is a generic prediction of our scenario that for profiles with 7 ~ 1 at r < r the gamma 
ray signal looks like DM decay from the inner region of the galaxy (cf. Eq. [T5"j) . whereas it 
mimics DM annihilation when looking away from the galactic center. 

Let us now discuss the dependence of the J-factor, Eq. [THl on A. In contrast to the usual 
annihilation case, J encodes not only the astrophysical DM profile, but it depends now also 
on the particle physics parameter A via p c g. This dependence is shown in Fig. [3] for the GC 
and GR regions observed by HESS, see Eqs. [H] and [TUl The figure shows that J gets reduced 
as soon as A becomes larger than the size x of the observed region, since then most of the 
<f) decay outside of the observed region. For the GC we have x ~ r & 'AQqc/tt w 15 pc, 
whereas for the GR we have x ~ r V A/A6 « 73 pc, in agreement with the values of A 
where the curves in Fig. [3] start to deviate from the J-factor for A = 0. In Tab. [I] we give 
J for the GC and GR for some values of A. In the case of GR the J factors also take into 
account that HESS subtracted the background by comparing to a region outside the center, 
as mentioned after Eq. [TTJ The quoted J factors are then a difference of the average fluxes 
from the two regions. For nonzero A this can lead to a further suppression of a factor of few. 
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III. APPLICATION TO RECENT COSMIC RAY DATA 



In this section we apply our LLP scenario to the recent data from cosmic ray experiments 
and discuss various other constraints. We describe our fit to the electron and positron data 
from PAMELA, ATIC, FERMI and HESS in Sec. ITTTAl and show that such a fit can be 
made consistent with the HESS GC and GR observations in Sec. IIIIB1 The bounds from 
neutrinos are discussed in Sec. IIII CI 



A. The electron-positron signals 

The electron flux from DM annihilations (which is equal to the positron flux) at r Q is 
calculated as 



dE e 4nb(E e 



M^jT M^IME.,*)), (23) 

with b(E e ) = E%/(GeV te), t e = 10 16 s, and the electron velocity v e ps c. dN e /dE is 
the injection spectrum for electrons (equal to the one for positrons) per <p decay, which we 
calculate by assuming that decays into muons. Then dN e / dE is calculated from the decay 
of the muons by using pythia-6.4. 19 [33( taking into account final state radiation. We 
provide analytic parameterizations of the injection spectra in Appendix |A] The diffusion 
length X D is given by X 2 D (E, E') = AK a T E (E & - 1 - E' 8 ~ l )/{1 - 5) with E, E' in GeV and 
throughout this work we assume the so-called MED propagation model from [341 ]. where 
K = 0.0112 kpc 2 /Myr and 5 = 0.70. The halo function I{Xd) is obtained as a series of 
Bessel- and Fourier transforms of Peff( r )/Peff( r o) m ^3]- We instead write down a 



partial differential equation for I(Xd) an d solve for it numerically, which speeds 
up the computation greatly. Details are relegated to appendix B. Note that the 
flux in Eq. [23] is a factor 2 larger than in the case of x annihilations directly into muons, 
since we obtain 2 0's for each XX annihilation, each of them giving a pair. 

We consider the measurement of the positron fraction $ e + / ($ e + + $ e - ) from PAMELA [l[ , 
where we use only the 9 data points above 6.8 GeV where the effect of solar modulation is 
expected to be small. Then we use data on the sum of electrons and positrons ($ e + + $ e -) 
from ATIC Q, FERMI jj, and HESS 0, Q (the effect of older data with much larger 
errors, e.g. PPB-BETS |35|, is expected to be small, therefore they are not included in the 
fits). Since ATIC and FERMI data are not consistent with each other we do not combine 
them but present results using only either of the two. For ATIC we use the combined data 
from ATIC 1, 2, and 4, see last reference in [3]. The two HESS measurements [f| and 0] 
overlap in the intermediate energy range. Therefore we use only the first 4 data points from 
0], while in the overlap region we take only the measurements of jfjj], which have smaller 
errors. The overall energy scale of FERMI and the two HESS data sets is varied within the 
quoted uncertainties. When we fit "PAMELA only" also the lowest three data points from 
ATIC are included in order to constrain the normalization of the background fluxes. For the 



astrophysical electron and positron background fluxes we use the parameterization from 36 



for the fluxes from Galprop [37j]. Following j38j, these fluxes are multiplied by C e ±E 
where in the fit we allow free normalization constants C e ± and assume a e ± = ± 0.05 (lcr), 
independently for e~ and e + . 
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FIG. 4: Allowed regions at 3ct for PAMELA (gray), PAMELA+HESS+FERMI (red), and 
PAMELA+HESS+FERMI (dark red) and the constraints from HESS photon observations of the galac- 
tic center (GC) and galactic ridge (GR) for A = 0, 1, 10, 100 kpc. The solid (dashed) curves show HESS 
photon constraints at 90% CL with (at 3cr without) including a power law background in the fit, see text 
for details. The regions above the curves are excluded. 



A [kpc] 


PAMELA+HESS+ATIC 
m x [TeV] (av) [chiV 1 ] x min 


PAMELA 

m x [TeV] ( 


+HESS+FERMI 

av) [cmV 1 ] Xmin 


PAMELA-only 
m x [TeV] (av) [cm 3 s _1 ] 


X 2 - 

Amin 





3.6 


1.2 x 10" 22 102.0 


3.5 


7.9 x 10~ 23 76.6 


0.40 


3.2 x 10~ 24 


7.3 


1 


3.6 


1.1 x 10~ 22 102.1 


3.5 


7.9 x 10~ 23 75.2 


0.35 


2.5 x 10~ 24 


7.1 


10 


3.5 


7.9 x 10~ 23 112.9 


3.6 


6.9 x 10" 23 82.9 


0.40 


2.5 x 10~ 24 


8.0 


100 


3.7 


6.3 x 10" 22 89.1 


3.7 


5.0 x 10~ 22 86.2 


0.32 


10 x 10" 24 


7.9 



TABLE II: Best fit values of m x and (av) and the corresponding Xj 2 nin values for electron-positron data 
for some representative values of the (j> decay length A. The number of degrees of freedom are 38, 43, 10 for 
PAMELA+HESS+ATIC, PAMELA+HESS+FERMI, PAMELA-only, respectively 



10 



I I I I III I I I I I IT T] 1 I I I III 1 I I I 1 1 Li 

A =10 kpc 

X = 

• PAMELA 



a 

o 



a 

o 

H 



o 
a. 




0.01 



10" 



10 10 
E [GeV] 



10 



I I I I I Mil 1 I I I Hill 1 I I I I Mil 1 I I I I ILL 



• PAMELA 

A =10 kpc 

X = 



a 

Q 

o 



o 

D. 




0.01 



10 io J 

£ [GeV] 



l(f 



a io 



S3 



10 



io J 



no ATIC 



I I I I I I I I I I — I I I I I 1 1 1 1 — I Mill 

• FERMI 
A ATIC 1-2-4 

HESS-08 

♦ HESS-09 




10 



10 

E [GeV] 



10 



a io z 



io 1 



n — i i i 1 1 1 1 i i — i i i 1 1 1 1 1 — i — tt" 

" 1 no FERMI 



FERMI 
ATIC 1-2-4 
HESS-08 
HESS-09 



" // X 2 =H2/38 
X 2 =102/38 

/ v I . I I 




J_l I l_L_ 



io z 



io J 

E [GeV] 



FIG. 5: PAMELA data on the positron fraction (left) and the electron-positron data (right) compared to 
the predicted spectra for A = (blue curves) and A = 10 kpc (black curves) at the best fit values given in 
Tab.lITl Solid curves correspond to signal + background, whereas with the dashed curves we show background 
and signal (right panels only) components separately. Upper panels are for PAMELA+FERMI+HESS, lower 
panels for PAMELA+ATIC+HESS. The green curves in the left plots show the spectrum at the best fit to 
only PAMELA data for A = 10 kpc. 



The results of our fit to electron data are shown as the shaded regions in Fig. H] in the 
plane of the DM mass m x and the x annihilation cross section (av) for four choices of the <p 
decay length A. These regions are defined by contours of A% 2 = 11.8 with respect to the \ 2 
minimum (3a for 2 dof). Since FERMI and ATIC data are inconsistent at about 3 a level, 
we perform separate fits where one of the two data sets is excluded. The corresponding best 
fit values are given in Tab. [Til an d the fit to the data is shown in Fig. [5j 

Irrespective of whether FERMI or ATIC data is used, the steep slope of the high-energy 
HESS data provides a strong constraint on the DM mass. This then leads to only a small 
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change for the best fit values of m x and (av) in the two cases. If we fit only PAMELA data, 
however, there is a degeneracy between m x and (av) , and as long as the peak in the electron 
spectrum is above the last PAMELA data point at 83 GeV a good fit is obtained. As shown 
in Fig. H] these results are basically independent of the decay length A, as long as this is 
not much larger than the distance from us to the center of the galaxy. The reason is that 
electrons and positrons are trapped in the turbulent galactic magnetic field and the observed 
signal is dominated by sources "near by" , within several kpc. Therefore, they are not very 
sensitive to the change of the source distribution from p 2 (r) to p 2 ff (r), which is mostly 
important close to the galactic center. Only if A becomes much larger than r « 8.5 kpc 
the total electron production close to us will be suppressed, which would require an increase 
in the annihilation cross section to maintain the signal. This effect is visible in Fig. H] and 
Tab. HTlfrom the results for A = 100 kpc. 



B. Gamma ray constraints from HESS GC and GR observations 



The photon fluxes predicted for given particle physics parameters m x , (av), A by Eq. [18] 
can be compared with observations. We use the gamma ray data from HESS observations 
of the GC and the GR Q. These data range from about 200 GeV to 20 TeV and are 
consistent with a power law spectrum oc E~ a with a ~ 2.3. In order to obtain bounds on DM 
parameters we adopt two different strategies. Most conservative bounds can be obtained by 
requiring that the signal predicted by DM must not exceed any data point of the observed 
flux. We obtain these bounds (denoted by "without background") by excluding points in 
the parameter space where the prediction exceeds x + 3a for any data point, where x is the 
observed flux and a its error bar. This leads to very conservative bounds, since it requires 
that in the signal region there is no astrophysical background, whereas in order to account 
for the observed flux at energies where DM does not contribute some astrophysical source 
has to be assumed. Therefore, we show also bounds by using a second method, called "with 
background". Here we fit the data with a power law background (with free normalization 
and power) + the signal from DM. For given A exclusion limits in the plane of m x and (av) 
at 90% CL are obtained by the contours with A% 2 = 4.6 with respect to the \ 2 minimum. 
The photon spectrum d/V 7 / dE^ per decay used in Eq. [TBI is calculated with pythia-6 .4.19 
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assuming that decays into fi + n~ + nj, see Appendix IA1 
The bounds from the HESS gamma ray observations are shown in Fig. H] together with 
the regions favored by electron-positron data. Clearly, for A = the electron data are 
inconsistent with the gamma ray constraints. The photon flux from the galactic center 
gets reduced for finite A, see Fig. EH and the bounds shift to larger values of (av) as A is 
increased. Fig. H] shows that for A ~ 10 kpc the regions favored by electron-positron data 
are consistent with the HESS gamma ray constraints, even for the 90% CL bounds with 
background. As one might expect the bounds are also consistent with purely decaying dark 
matter 3ll. l39l . |40| (and even in the case where decaying dark matter is also allowed to 
annihilate [411]). 

Apart from these high energy gamma ray observations additional information can be 
obtained from less energetic photons. At energies around 10 GeV the photon flux comes 
mainly from inverse Compton scattering on CMB photons, star-light and dust-light, both 
for annihilating 42, 43, 44|, 45] and decaying DM 46, 47, 48 1. This leads to a diffuse gamma 
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ray signal, potentially observable by FERMI in the near future. The recent diffuse gamma 
flux measurement by FERMI in the 10° — 20° band above the galactic plane leads to a 
bound an order of magnitude weaker than needed to probe the m x — (av) region preferred 
by electron/positron data even for the case of A = (22)]. Since the effective density is 
smaller in the inner galactic region for A > 0, cf. Fig. [TJ this bound is expected to become 
even weaker for the LLP scenario. Looking away from the inner galactic region the signal 
can, however, be stronger than for the standard annihilating DM. To asses the impact of 
FERMI a more detailed analysis is called for. 

In addition to constraints from gamma rays, there are also strong constraints due to 



synchrotron radiation of radio waves. These bounds were extensively studied in ly, [18j, 



49l . |5(| where it was found that annihilation explanations of the positrons seem to be ruled 



out, unless the DM profiles are made effectively less steep as in decays or as in LLP. Thus 
it seems very plausible that for sufficiently large A LLP can satisfy the bounds, though 
calculations along the lines of ly, |l_8j are in order. Constraints usin g po tential gamma ray 



signals from dwarf galaxies, on the other hand, are less powerful [51|, 152 
galaxy clusters have been discussed in [53 . 



Constraints from 



C. Neutrinos from the galactic center 



Super-Kamiokande (SK) provides an upper limit on the upward going muon flux from 



various extra-terrestrial sources 54 



sec 



55 for similar recent results from Ice Cube. Since 



the LLP scenario predicts 8 neutrinos for each DM annihilation (from the decay of the four 
muons from XX ~ ¥ 00 ~ * 2/^ + 2/x~) the bound on upward going muons from the galactic 
center is potentially relevant. The neutrino induced muon flux can be calculated as, see 
e.g. 

dxj%{E v ) 



/ dE »7W dE,R,(E,)J2n a J2 

JE ihl CSA/ J Ethl . a=p n x=v 9 



dE„ 



Here, 



In 



(24) 



(25) 



is the range of a muon with energy until its energy drops below £7thr> f° r which we take the 
SK analysis threshold of 1.6 GeV, with a M = 2x 1(T 3 GeVcm 2 g" 1 , (3^ = 4.2 x 1CT 6 cm 2 g _1 , 
and p is the density of the material in gem -3 . Further, n a « r a p/m p are the number 
densities of neutrons and protons with r p « 5/9, r n « 4/9, and for the detection cross 
section we use 

do%(E v ) 



dE„ 



2m p G 2 F 

71 



El 



(26) 



with Ay = 0.25, 0.15, By = 0.06, 0.04, and A n / = By, B^ p = Ay [57 . 

dQ v /dE v is the flux of muon neutrinos arriving at the earth within a solid angle AQ. In 
our case the flux is equal for neutrinos and anti-neutrinos, and it is given by 



d$„ = (av) r Q p 2 Q 
dE v 47r ml 



dN u , 
dE v 



dN Uf 
dE v 



JAtt 



(27) 



13 




The oscillation probabilities in terms of the lepton mixing matrix elements U a i are P Ve -^ v = 

Etil^eil 2 !^! 2 « °- 21 > = £-=i ~ °- 395 - dN UeM /dE u is the spectrum of 

electron (muon) neutrinos per <p decay (and integrates to 1). A parameterization for it is 
given in Appendix |A] The J-factor in Eq. [27] is defined in the same way as for photons, see 
Eq. [TjH and therefore the neutrino flux will be reduced with increasing the (f) decay length A 
similar to the case of photons. 

SK provides an upper bound on the muon flux from a cone around the galactic center 
as a function of the cone half opening angle up to 30°. In Fig. [6] we compare the SK upper 
bound to the predicted muon fluxes for fixed m x and (av) as a function of A. We conclude 
that A = is close to the present bound, while for A ^ 5 kpc the SK bound is significantly 
relaxed. For SK neutrino bounds for decaying DM see 58|, |59|, [60 . 

Let us also mention that in the LLP scenario we do not expect any observable neutrino flux 
from DM annihilations in the sun (or in the earth), since <ft is to a very good approximation 
stable at the scale of the solar system. These neutrino fluxes will be exponentially suppressed 
by the ratio of the sun-earth distance (or the earth radius) to A. 



IV. CAN x BE A THERMAL RELIC? 

Let us now address the issue of whether or not \ can be a thermal relic. The LLP 
scenario, where x annihilates first to 2 0's and these then after some time decay to SM 
particles, faces similar challenges as direct x' s annihilation into SM particles. 

First of all, for x to be a thermal relic, the required \ annihilation cross section at the time 
of freeze-out in the early universe is several orders of magnitudes too small to explain the 
PAMELA/ATIC/FERMI cosmic ray anomaly. For x with TeV mass, the cross section at the 
freeze-out should be about (av )fo ~ 3 x 10~ 26 cm 3 s~ 1 , while as we have seen in section UlI A] 
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the annihilation cross section that explains the cosmic ray anomaly, is much larger, (<jv)pa ~ 

1 One intriguing possibility which can explain the mismatch is that there exists 

"1 2 



10 22 cm 3 s 



an attractive long range force between the dark matter particles [6JJ • This leads to an 
enhancement of the cross section at small velocities by a factor of l/v. With a typical 
velocity of x's i n the galactic halo v/c ~ 10 -3 this so called "Sommerfeld enhancement" is 
roughly of the right size. In the LLP scenario for the case of a non-relativistic 0, there is 
a final phase space suppression of the cross section for x's annihilating in the galactic halo 
compared to the early universe. However, larger Sommerfeld enhancements are attainable, 
if there exists a bound state very close to threshold [62|, [63 . 

The Sommerfeld enhancement solution has several potential phenomenological problems. 
For instance, if the Sommerfeld enhancement worked to arbitrarily small velocities, then 
annihilation in proto-halos could lead to a too large contribution to the diffuse gamma ray 
background [64J. For massive enou gh attractive force carriers this bound does not apply 
(such as GeV mass force carriers of 24j) as the enhancement saturates. More importantly, 
the highly energetic leptons and photons originating from x annihilations at T < 10 keV in 
the early universe could lead to photo-dissociation of light elements destroying the successful 
predictions of standard big bang nucleosynthesis (BBN) j65[. The resulting bound on the 



annihilation cross section for XX ~> A* 4 A* * s ( av )FO < 2.0 x 10 



-23 cm 3 s 



(m x /lTeV)" 



This bound is already somewhat smaller than needed for the explanation of the cosmic ray 
anomaly, see FigHl The bound itself was obtained assuming time independent cross section, 
so a more detailed study in the framework of models giving Sommerfeld enhancement (and 
in the LLP scenario) may be warranted. However, since in our case cf) is very long lived 
compared to BBN time scales, one may expect that within the LLP scenario there is no 
threat for BBN from rapid x annihilations, while there are constraints from the late decays 
of </>'s, see below. 

There are other ways to "boost" the annihilation cross section in the galactic halo with- 
out running into these phenomenological problems. For instance the annihilation could go 



through a resonance with mass of order 2m x 66|, |67|. This could arise naturally in models of 



extra dimensions with linearly spaced KK modes. Another possibility would be a kination 



model 68|, [69], where the expansion rate at the time of decoupling is increased due to a 



rolling scalar field, leading to a reduced relic density. Yet another possibility is that x is 



a product of a decay of a different meta-stable thermal relic 70j. Clearly, no boost factor 



is needed for decaying DM, but the large decay time leads to interesting model-building 



implications [53, |7l|, |72j, [73 



In the LLP scenario, where x is a thermal relic and <fi a meta-stable thermal relic, there 
are other constraints coming from the late decay of the particle. Phenomenologically, the 
LLP scenario is interesting if travels a distance A ^ 10 kpc, because then it has a significant 
impact on the expected photon flux from the galactic center as discussed in section IIII Bl 
This means that the <b life time is 



A 



10 12 : 



c/?7 /?7 \ 10 kpc 



A 



(28) 



Such late decaying relics can lead to modifications of the light element abundances [74], [75 



Long range means here that the mass of the force carrier is smaller than Mv. 
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For TeV masses and lifetimes in this range, the particle must have a relic density which 

For 



is three orders of magnitude smaller than the one for \i see figures 5 and 6 of [75[ 
t > 10 13 s there is an even stronger bound on the abundance due to constraints on 



the diffuse gamma ray background [74|, |76| (see figure 7 of [76|, where the bounds were 



derived for the case of radiative decays, while the bounds for leptonic decays are expected 
to be slightly weakened). The ratio a^ nn /<70 nn must therefore be of order 10~ 3 or smaller. 
Consequently, if the couplings to quarks are not excessively suppressed, pair production 
should be observable at LHC We leave it as an open exercise in model building to embed 
the above hierarchy of cross sections in a concrete model. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we introduced a scenario for cosmic ray production due to dark matter 
(DM) which is an alternative to annihilations or decays 3 . The scenario softens the angular 
distribution of the gamma ray spectrum by having the DM annihilate into long lived particles 
(LLP), that subsequently decay to SM particles. We have demonstrated within the context of 
the PAMELA/ATIC/FERMI positron and electron cosmic ray anomaly, that while standard 
DM annihilation scenarios seem to be in conflict with the gamma ray data due to large 
required cross sections {{<jv) ~ 10 -23 cm 3 s _1 ), the prolongation of the final state lifetimes in 
the LLP scenario can accommodate the data. It should be emphasized that the gamma ray 
bounds are sensitive to the dark matter profiles. For less cuspy profiles the bounds would 
be weakened. Here we have only considered the NFW profile which grows as 1/r for small 



r. Recent simulations [77J, [78( that include baryons seem to indicate an even faster rise for 
small r which would only strengthen the gamma ray bounds. 

An attractive feature of DM annihilations is the direct relation between the annihilation 
cross section for indirect DM detection and the production cross section at collider exper- 
iments. In the LLP scenario this connection is lost to some extent, since the annihilation 
cross section of x is n °t directly related to the production cross section of 0. However, 
assuming that x an d are thermal relics, bounds from light element abundances and diffuse 
gamma rays on decaying relics imply that the annihilation cross section of has to be some 
orders of magnitude larger than the one of x 111 order to suppress its relic abundance, see 
discussion in section IIVI Under the assumption of a non-negligible coupling to hadrons this 
implies good prospects for production at LHC (if kinematically accessible). 

It is interesting to ask how one would differentiate the LLP scenario from DM decay. The 
FERMI satellite should be able to accurately measure the angular distribution of the gamma 
ray flux. This distribution can be used to differentiate between decays and the LLP scenario 
introduced in this paper. Figure [7] shows the gamma ray flux at fixed galactic latitude b = 15 
degrees as a function of the galactic longitude for the decay and the annihilation profiles 
together with the one for the LLP scenario with A = 10 kpc. We see that while the long 
lifetime reduces the slope compared to annihilation it still cuts off faster than the decay 
profile. An analysis of the angular distribution for decaying dark matter was performed in 
[39( and it would be interesting to extend this analysis to our scenario. 



A scenario where the dark matter both anhinilates and decays was propsosed in 41 1 . 
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FIG. 7: Angular distribution of the gamma ray flux for the LLP scenario with A = 10 kpc in comparison 
to the ones for DM decay and annihilation. We show the gamma ray flux for fixed galactic latitude 6=15° 
as a function of the longitude I, normalized to its value at 30°. 

Another distinguishing feature of LLP will be the smearing out of clumps. If the only 
boost factor were due to clumps then this would be a striking signal for dark matter. Typi- 
cally the over-densities lie towards the edges of galaxies where they are not tidally disrupted. 
The long lived particle scenario would smooth out this striking signal. 

Here we have explored only the case of non-relativistic intermediate states, which implies 
a near degeneracy between the DM and the LLP. However, the softening of the angular dis- 
tribution demonstrated here, which allows one to avoid constraints from gamma ray bounds, 
should be active in the case of relativistic annihilation products as well. An investigation of 
this case is warranted. 

Finally, let us speculate on the possible origin of the energy scale corresponding to the 
typical life times of the LLP in our scenario, which is set by our distance to the galactic 
center of order 10 kpc, see Eq. [2SJ If we assume that the decay of the LLP <p is governed 
by a dimension six operator suppressed by a high scale A one has roughly r ~ 167rA 4 /m0. 
Then we find from Eq. for the scale A 



masses, since v 2 /A ~ 1 eV for v ~ 100 GeV and A ~ 10 13 GeV. 4 This is particularly 
intriguing in light of the leptonic nature of the currently observed cosmic ray features and 
might indicate some common high scale physics responsible for light neutrino masses and 
the LLP decay in our scenario. 



4 We mention that a similar argument for a dimension five operator would point to a scale A above the 




(29) 



We observe that the scale 10 13 GeV corresponds roughly to the seesaw scale for neutrino 



Planck scale. 
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APPENDIX A: INJECTION SPECTRA 



Here we provide injection spectra for photons, electrons, and neutrinos resulting from the 
decay <fi — > fi + n~ in 0's rest frame. 5 Since we work in the limit where 0's are non-relativistic, 
the provided injection spectra can be directly used in the numerical analysis in section IIII1 
The spectra are obtained by fitting analytic functions to a Monte Carlo sample generated 
with pythia-6.4. 19 33j. In the following we define x = lEjm^ and < x < 1. 



Photon injection spectra. The photon spectrum for the decay (or similarly annihi- 
lation) of a heavy particle with mass into charged particles is well described by the 
Weizsacker- Williams form (see e.g. 171. Il8l|) 



dN 1 
dK 



EyTl 



X) 



1] In 



ml(l - x) 



mf 



(Al) 



To be specific we have assumed here decay into a lepton pair /. Though in our work we 
consider only the muon final state we list here for completeness also the photon spectrum 
for electron and tau final states. The results from Pythia for muon and electron final states 
are very close to Eq. IAU In the range for E 1 from 10 GeV to 500 GeV for 1 TeV decaying 
DM the Pythia spectrum is 0.8 to 0.6 of the Weizsacker- Williams form. We parameterize 
the spectrum by 



dK 



A 



Or 



E y TT 





V4(l - x)~ 


(in 







+ B) [{l- x f + C}. 



(A2) 



The coefficients in the above formula depend on the value of m^. We parameterize this 
dependence by 

K{m * ) = Ko (mkT ' (A3) 

where K = A, B, C. Fitting Eq. IA2I to the spectrum generated with Pythia we find for 
mj, e [10 2 , 10 4 ] GeV: 



e + e final state: 
final state: 



(A Q , Bo, Co) = (1.49197, -11.0576, 0.688768) 
(5 A , S B , 5c) = (0.050141, 0.235939, -0.00623121) ' 

(Aq, Bo, Co) = (1.36693, -2.94774, 0.517128) 
(5 A , 5 B , 8c) = (-0.0365572, 0.412598, -0.0137187) 



(A4) 



5 The spectra apply also for the case of standard annihilations \X — * M + M by replacing — > 2m x . 
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The photon spectrum is then described to better than 13% (18% for /i + /i~) everywhere in 
the fitted range (mostly to better than 5%). This was obtained from 6 values of m^, for 
each of the values with 10 7 simulated decay events. 

For the decay into t + t - we have produced 10 7 events in Pythia. Since the taus are heavy 
and can decay, producing hadrons and leptons, and these further decay or bremsstrahl 
photons, the Weizsacker- Williams form Eq. IA1I is no longer a good description of the final 
photon spectrum. A form that gives a good description is 

dN y a cm r , 2 . B B 1 

M, = A °E^ eXp[ - (AlX + A2X)] + lGeV + IGeV* ' (A5) 



with Ai, Bi depending as in Eq. IA3I on with 

(Aq, A\, A° 2 , B\, Bl) = (283.066, -0.981786, 5.80398, -0.000155047, 0.000145457) 
(5 Ao , S Al ,5 A2 , S Bl ,S B2 ) = (0.00118045, -0.0141779, -0.00573809, -0.966239, -0.964131) 

(A6) 

The error on the predicted spectra is below 15% (in the first bin it is 30%). 

Electron injection spectra. Consider a two-body decay <p —* and neglect the 

photon emission. Then the electron spectrum can be obtained by boosting the isotropic 
spectrum from the muon rest frame into the rest frame of (p. In our work we do include final 
state radiation and therefore, we prefer to use the electron spectrum generated with Pythia. 
We use the function 

M J^ n v , = ^o— exp[-(A 1 x + A 2 x 2 )] + B + B x x , (A7) 
d(E e /l GeV) 7r 



with Ai, Bi depending as in Eq. [A3] on m^. Fitting this form to the generated spectrum we 
find 

(A° , A\, A° 2 , B°, B°) = (-0.296635, 2.65121, 14.8445, 0.0042505, -0.00427157) , , 
(6 Aoi 6 Al ,5 A2 ,6 Boi 6 B J = (-1.01424, 0.017198, -0.0107585, -0.999819, -0.999819) ' ^ ' 

The error on the predicted spectra is below 7%. 

Neutrino injection spectra. In each decay <p —* [i + fi~ we have two neutrinos — and 
v e — and two antineutrinos with the same spectra. It turns out to be convenient to fit for 
the full neutrinos spectrum (y e + v^) and for the electron neutrino spectrum separately. The 
muon neutrino spectrum is then the difference between the total neutrino spectrum and the 
electron neutrino spectrum. For the total v e + neutrino spectrum (normalized to 2) we 
use the parameterization 

_^_.^ exph( ^ + ^ )]+B , (A9) 

with Ai, B depending as in Eq. [A3] on and 

(A° , Al, A% B°) = (3.73892, 0.0128514, 2.0015, -0.00117002) 
(5 Ao , 5 Ai ,5 m ,5b) = (-1.0019, -0.283851, 0.00350415, -1.01054) 



(A10) 
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with the error on the predicted spectra below 4% (but mostly below 2%). And for the 
electron neutrino spectrum we use 

diEjrkev) = Ao ^ exp + Ax2 + A ^^^)] + B > ( An ) 

with again A+, B depending as in Eq. IA3l on m^, and 

(A° , A?, A° 2 , B°, A° 1/2 ) = (14.4735, 0.90821, 2.90886, -0.000653652, 1.95634) 
{6 Ao , S Al ,S A2 ,S B , 5 Al/2 ) = (-0.994217, -0.00194938, 0.00048921, -0.985772, -0.000395335) 

(A12) 

The v e spectrum is described to better than 5%. 



APPENDIX B: DIFFERENTIAL EQUATION FOR I(X D ) 

The propagation of positrons in the galaxy is described by a transport equation (see, e.g. 

H) 

d t V - V(K(x,£)W) -d E (b(E)V) = q(x,E), (Bl) 

with ^(x, E) the positron number density per energy interval and q(x, E) the positron source 
term. The loss coefficient due to Compton scattering is b(E) = E € 2 /te where e = E/E 
with Eq = 1 GeV and te = 10 16 s, while the transport through the turbulent galactic 
magnetic field is approximated by a spatially constant diffusion coefficient K(x, E) = Kqc s . 
For annihilating dark matter the positron injection source term is 

q (x,E) = J^) 2 f(e), (B2) 

while for LLP scenario p e fj instead of p should be used (for decaying dark matter, on the 
other hand, the dependence on p is linear). The positrons flux at radial distance from 
galactic center r and distance z from the galactic plane is then given by 

POO 

*(r,z,e) = K^j de s f(e s )I(X D ,r,z), (B3) 

where the diffusion length is A^> = y/4K Q f, with f = t — t s and t = r s (e' 5_1 /(l — 5)). The 
function I(\£),r,z) depends only on astrophysics, i.e. only on the DM halo profile. The 
function /(A^) in Eq. [23] is then equal to I(Xi),r, z) at the position of the solar system 

I(X D ) =/(A D ,r ,O). (B4) 

In [34| the function J(Ai),r, z) is given in terms of an infinite sum over Bessel and cosine 
functions, however this expansion converges very slowly. 

Alternatively, using the ansatz Eq. IB3l in Eq. IB 1L one can write down a partial differential 
equation for / 

V 2 /(A D , r, z) - ^d x J(X D , r, z) = 0. (B5) 
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The boundary conditions 




1(0, r, z) 




(B6) 



0. 



follow from a requirement that ^(r, z, e) vanishes on the boundaries of the galactic disk with 
radius R ga \ and thickness 2L (the first boundary condition in Eq. IB6l has to be multiplied by 
a smooth step function near the galactic disk boundaries to make it consistent with the next 
two boundary conditions). The above partial differential equation can be solved numerically 
easily using standard methods, for instance the method of lines. The solution is obtained on 
a personal computer in less then a minute for the whole galaxy, compared to several hours 
for each point, if the expansion in Bessel and cosine functions is used. With appropriately 
modified boundary conditions the above partial differential equation can be used for LLP 
scenario as well as for decaying dark matter. 



[1] O. Adriani et al. larXiv:0810.4995l [astro-ph]. 

[2] J. Chang et al, Nature 456 (2008) 362; preliminary results from the ATIC-4 flight were shown 
by J. Isbert at the workshop "Tango in Paris", May 4-6, 2009, Institut d'Astrophysique de 
Paris, France. 

[3] A. A. Abdo et al. [The Fermi LAT Collaboration] , larfOviOgOLQQ^ [astro-ph. HE]. 

[4] A. Buffmgton, C. D. Orth and G. F. Smoot, Astrophys. J. 199, 669 (1975); D Mueller 

and K.K Tang, Astrophys. J. 312, 183 (1987); R. L. Golden et al, Astrophys. J. 436, 

769 (1994); S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482, L191 (1997) 

|arXiv:astro-ph/9703T92] . 
[5] S. Torrii et al, Astrophys. J. 559, 973 (2001); J. Nishimura et al, Advances in Space Research 

19, 767 (1997); J. Nishimura et al. Astrophys. J. 238, 394 (1980). 
[6] F. Aharonian et al. [H.E.S.S. Collaboration], Phys. Rev. Lett. 101, 261104 (2008) 

|arXiv:0811.3894l [astro-ph]] . 
[7] H. E. S. Aharonian. rarXiv:0905.0T05l [astro-ph.HE]. 
[8] J. H. Taylor and D. R. Stinebring, ARA&A 24, 285 (1986). 

[9] D. C. Backer, Pulsars (Galactic and Extragalactic Ra- dio Astronomy, 1988), pp. 480521. 



[10] D. Hooper, P. Blasi and P. D. Serpico, JCAP 0901, 025 (2009) |arXiv:0810. 15271 [astro-ph]]. 

[11] H. Yuksel, M. D. Kistler and T. Stanev, arXiv:0810.2"784l [astro-ph]. 

[12] S. Profumo. arXiv:0812.4457l [astro-ph]. 

[13] N. J. Shaviv, E. Nakar and T. Piran. larXTv:0902.0376l [astro-ph.HE]. 

[14] D. Malyshev, I. Cholis and J. Gelfand. larXiv:0903. 13101 [astro-ph.HE]. 

[15] J. Hall and D. Hooper. la7Xiv:0811.3362 [astro-ph]. 

[16] G. Bertone, M. Cirelli, A. Strumia and M. Taoso. larXiv:0811.3744l [astro-ph]. 

[17] N. F. Bell and T. D. Jacques. [arXiv:0811.082i1 [astro-ph], 

[18] L. Bergstrom, G. Bertone, T. Bringmann, J. Edsjo and M. Taoso, arXiv:0812.3895 [astro-ph]. 



(2008). 



21 



[19] P. Meade, M. Papucci and T. Volanskv. la"rXiv:0901.2 925 [hep-ph]. 

[20] F. Aharonian et al. [H.E.S.S. Collaboration], Phys. Rev. Lett. 97, 221102 (2006) [Erratum- 

ibid. 97, 249901 (2006)] |arXiv:astro-ph/0610509| . 
[21] F. Aharonian et al. [H.E.S.S. Collaboration], Nature 439, 695 (2006) |arXiv:astro-ph/060302i] . 
[22] P. Meade, M. Papucci, A. Strumia and T. Volanskv. la"rXiv:0905.0480l [hep-ph]. 
[23] M. Pospelov and A. Ritz, Phys. Lett. B 671, 391 (2009) [arXiv:0810.1502l [hep-ph]]. 
[24] N. Arkani-Hamed, D. P. Finkbeiner, T. Slatyer and N. Weiner, Phys. Rev. D 79, 015014 

(2009) [arXiv:0810.0713l [hep-ph]]. 
[25] I. Cholis, D. P. Finkbeiner, L. Goodenough and N. Weiner. larXiv:0810.5344l [astro-ph]. 
[26] Y. Nomura and J. Thaler. larXiv:0810.5397l [hep-ph]. 

[27] I. Cholis, G. Dobler, D. P. Finkbeiner, L. Goodenough and N. Weiner. [arXiv:0811.364"T1 [astro- 

Ph]. 

[28] P. J. Fox and E. Poppitz. larXiv:081 1.03991 [hep-ph]. 

[29] F. Chen, J. M. Cline and A. R. Frev. iarXiv:0901 .43271 [hep-ph]. 

[30] J. Mardon, Y. Nomura, D. Stolarski and J. Thaler. larXiv:0901 .29261 [hep-ph]. 

[31] E. Nardi, F. Sannino and A. Strumia, JCAP 0901 (2009) 043 |arXiv:081 1.41531 [hep-ph]]. 

[32] J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 462, 563 (1996) 

|arXiv:astro-ph/9508025] . 
[33] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006) |arXiv:hep-ph/0603175| . 
[34] T. Delahaye, R. Lineros, F. Donato, N. Fornengo and P. Salati, Phys. Rev. D 77 (2008) 063527 

|arXiv:0712.23T2l [astro-ph]]. 
[35] S. Torii et al. larXiv:0809.0760l [astro-ph], 

[36] E. A. Baltz and J. Edsjo, Phys. Rev. D 59, 023511 (1999) |arXiv:astro-ph/9808243| . 
[37] I. V. Moskalenko and A. W. Strong, Astrophys. J. 493 (1998) 694 |arXiv:astro-ph/9710124~ . 
[38] M. Cirelli, R. Franceschini and A. Strumia, Nucl. Phys. B 800, 204 (2008) |arXiv:0802.3378l 
[hep-ph]]. 

[39] G. Bertone, W. Buchmuller, L. Covi and A. Ibarra, JCAP 0711, 003 (2007) [ arXiv:0709.2299l 
[astro-ph]]. 

[40] A. Ibarra and D. Tran, larXiv:0811.T555l [hep-ph]. 

[41] K. Cheung, P. Y. Tseng and T. C. Yuan. larXiv:0902.4035l [hep-ph], 

[42] C. R. Chen, M. M. Nojiri, S. C. Park, J. Shu and M. Takeuchi. [arXiv:0903.197ll [hep-ph]. 
[43] E. Borriello, A. Cuoco and G. Miele. larXiv:0903.1852l [astro-ph.GA], 
[44] M. Lattanzi and J. I. Silk. larXiv:0812.0360 [astro-ph]. 
[45] R. Harnik and G. D. Kribs. larXiv:0810.5557 [hep-ph]. 
[46] K. Ishiwata, S. Matsumoto and T. Moroi. larXiv: 0903 .02421 [hep-ph], 
[47] S. Shirai, F. Takahashi and T. T. Yanadda. larXiv:0902.4770l [hep-ph]. 
[48] P. Grajek, G. Kane, D. Phalen, A. Pierce and S. Watson. larXiv:0812.4555l [hep-ph]. 
[49] J. Zhang, X. J. Bi, J. Liu, S. M. Liu, P. f. Yin, Q. Yuan and S. H. Zhu, arXiv:0812.0522l 
[astro-ph] . 

[50] K. Ishiwata, S. Matsumoto and T. Moroi, Phys. Rev. D 79, 043527 (2009) |arXiv:081 1.44921 
[astro-ph]]. 

[51] R. Essig, N. Sehgal and L. E. Strigari. larXiv:0902.4750l [hep-ph], 
[52] L. Pieri, M. Lattanzi and J. Silk. larXiv:0902.4330l [astro-ph.HE]. 
[53] Q. Yuan, X. Bi, F. Huang and X. Chen. larXiv:0902.4294l [astro-ph.CO]. 



22 



[54] S. Desai et al. [Super-Kamiokande Collaboration], Phys. Rev. D 70, 083523 (2004) [Erratum- 
ibid. D 70, 109901 (2004)] | arXiv:hep-ex/0404025| . 
R. Abbasi, larXiv: 0902 .24601 [astro-ph.CO]. 
C. Delaunay, P. J. Fox and G. Perez, arXiv:0812.333"Tl [hep-ph]. 

V. Barger, W. Y. Keung, G. Shaughnessy and A. Tregre, Phys. Rev. D 76, 095008 (2007) 
|arXiv:0708. 13251 [hep-ph]]. 

J. Hisano, M. Kawasaki, K. Kohri and K. Nakayama, arXiv:0812.0219 [hep-ph]. 
A. Arvanitaki, S. Dimopoulos, S. Dubovsky, P. W. Graham, R. Harnik and S. Rajendran, 
larXiv:Q812.2075l [hep-ph]. 

J. Liu, P. f. Yin and S. h. Zhu. larXiv:0812.0964 [astro-ph]. 

J. Hisano, S. Matsumoto, M. Nagai, O. Saito and M. Senami, Phys. Lett. B 646, 34 (2007) 
[arXiv:hep-ph/0610249] . 

J. D. March-Russell and S. M. West, arXiv:0812.0559l [astro-ph]. 

M. Pospelov, A. Ritz and M. B. Voloshin, Phys. Lett. B 662, 53 (2008) |arXiv:071 1.4 866 
[hep-ph]]. 

M. Kamionkowski and S. Profumo, Phys. Rev. Lett. 101, 261301 (2008) [arXiv:0810.3233l 
[astro-ph]]. 

J. Hisano, M. Kawasaki, K. Kohri, T. Moroi and K. Nakayama, arXiv:0901.3582 [hep-ph]. 
M. Ibe, H. Murayama and T. T. Yanagida, larXiv:0812.0072l [hep-ph]. 
W. L. Guo and Y. L. Wu. larXiv:0901.1450l [hep-ph]. 
P. Salati, Phys. Lett. B 571, 121 (2003) [ arXiv:astro-ph/0207396l . 

F. Rosati, Phys. Lett. B 570, 5 (2003) [arXiv:hep-ph/0302159| . 
M. Fairbairn and J. Zupan, arXiv:0810.4147 [hep-ph]. 
M. Pospelov and M. Trott, larXiv:0812.0432l [hep-ph]. 

P. f. Yin, Q. Yuan, J. Liu, J. Zhang, X. j. Bi and S. h. Zhu, Phys. Rev. D 79, 023512 (2009) 
|arXiv:0811.0T76l [hep-ph]]. 

C. R. Chen, M. M. Nojiri, F. Takahashi and T. T. Yanagida. larXiv:081 1.33571 [astro-ph]. 
J. R. Ellis, G. B. Gelmini, J. L. Lopez, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 373, 
399 (1992). 

M. Kawasaki and T. Moroi, Astrophys. J. 452, 506 (1995) [arXiv:astro-ph/9412055] . 

G. D. Kribs and I. Z. Rothstein, Phys. Rev. D 55, 4435 (1997) [Erratum-ibid. D 56, 1822 
(1997)] [arXiv:hep-ph/9610468] . 

O. Y. Gnedin, A. V. Kravtsov, A. A. Klypin and D. Nagai, Astrophys. J. 616, 16 (2004) 
|arXiv:astro-ph/0406247] . 

S. Kazantzidis et al, Astrophys. J. 623, L67 (2005) [arXiv:astro-ph/0407407] . 



23 



